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FOREWORD 


This, Volume II of the report, "Spherical Roller Bearing 
Analysis," details information required to use the design and 
analysis computer program SPHERBEAN. All efforts involved in 
the generation of the code were accomplished under contract 
NAS3-20824 issued by the NASA-Lewis Research Center of Cleveland, 
Ohio, under the administration of the Structures and Mechanical 
Technologies Division. Funding was provided by the Product 
Assurance Office of the Army Aviation Research and Development 
Command, St. Louis, Missouri. The technical monitor was 
Mr. H. H. Coe. The work was performed at SKF Industries, Inc., 
King of Prussia, Pennsylvania, during the period June, 1978 
through December 1980. 

Technical project leadership was executed by Mr. R. J. 
Kleckner, with contributions from: Dr. J. Pirvics, Messrs. 

G. Dyba, T. Deromedi and Ms. M. M. Dinon. 


SKF TECHNOLOGY SERVICES 

SKI- INOUSIPItS. INC 


TABLE OF CONTENTS 


SECTION 

1.0 

2.0 

3.0 

3.1 

3.2 

3.3 

3.4 

3.5 

3.6 

4.0 

4.1 

4.2 

5.0 

6.0 

APPENDIX A 
APPENDIX B 
APPENDIX C 
APPENDIX D 


TITLE PAGE 

INTRODUCTION 1 

INPUT DATA 3 

SPECIAL INPUT DATA 28 


OPTION 1: COMPUTATION OF BEARING 

OPERATING DIAMETRAL CLEARANCE . 29 

OPTION 2: USER SPECIFIED MATERIAL 


PROPERTIES 30 

OPTION 3: USER SPECIFIED SLICE WIDTHS ... 31 

OPTION 4: USER INPUT SYMMETRIC ROLLER 

GEOMETRY 32 

OPTION 5: USER SPECIFIED, COMPLETELY 

VARIABLE, ROLLER GEOMETRY 33 

OPTION 6: TEMPERATURE CALCULATIONS 34 

OUTPUT DATA 4 5 

EXAMPLE 1: SIMULATION OF STEADY STATE 

BEARING PERFORMANCE 46 

EXAMPLE 2: PLANETARY BEARING ANALYSIS ... 53 

PROGRAM LIMITATIONS AND GUIDES TO 

PROGRAM USE 54 

LIST OF REFERENCES 57 

INPUT DATA AND FORMATS A1 - A6 8 7 

HEAT TRANSFER COMPUTATION NOTES 106 

SPHERBEAN OUTPUT, EXAMPLE PROBLEM 1 114 

SPHERBEAN OUTPUT, EXAMPLE PROBLEM 2 139 




SKF TECHNOLOGY SERVICES 

INl JUSTICES. B'-j 



LIST OF TABLES 


NO. TITLE I’AOl: 

1. ORGANIZATION OF INPUT DATA CATEGORIES 58 

2. DEFAULT VALUES FOR USER UNSPECIFIED VARIABLES ... 59 

3. PROPERTIES OF FOUR LUBRICANTS 61 

4. ORGANIZATION OF SPECIAL INPUT DATA 6 2 

5. SPECIFICATIONS OF SPHERICAL ROLLER BEARING 

USED IN EXAMPLE 1 6 3 

6. DICTIONARY OF NODES USED IN SYSTEM MODEL, 

EXAMPLE 1 64 

7. TWO DIAGNOSTIC MESSAGES GENERATED IN THE 

NEWTON- RAPHSON EQUATION SOLVER 65 


l- II: Of (NGlOGY SLRVlCES 


i i 


LIST OF FIGURES 


NO. 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 
1 5 
10 


17 


18 

19 

20 
21 


TITLE PAGE 

SPHERICAL ROLLER BEARING 06 

USE OF FREE FORMAT TO SPECIFY INPUT DATA 67 

USER INPUT ROLLER GEOMETRY 08 

BEARING COORDINATE FRAME 69 

USER INPUT BEARING CLEARANCES 70 

USER INPUT RING DATA 71 

PLANET BEARING 72 

PLANETARY TRANSMISSION 73 

HIGH SPEED PLANET GEAR LOADING 74 

PLANET GEAR OUTER RING CROSS-SECTION 75 

IDEALIZED GEAR TOOTH LOADS 76 

TOOTH LOADS ON PLANET GEAR (NMESH =3) 77 

CAGE GEOMETRY 7 8 

SPLIT CAGE INPUT DATA 7 9 

TEMPERATURE NODE IDENTIFICATION SCHEME 80 

OPTIONAL SYMMETRIC ROLLER GEOMETRY INPUT DATA. 81 

OPTIONAL ROLLER GEOMETRY DATA 82 

SPHERICAL ROLLER BEARING TEST RIG 83 

SYSTEM MODEL USED IN EXAMPLE 1 (METAL AND 

AIR NODES) 8 4 

SYSTEM MODEL USED IN EXAMPLE 1 (LUBRICANT 

SYSTEM NODES) 85 

ROLLER COORDINATE FRAMES USED IN PROGRAM 

OUTPUT 86 


m ■mu'n jjihbjii jj.bi hjmp 






AT81D007 


SPHERICAL ROLLER BEARING ANALYSIS 
1.0 INTRODUCTION 

SPHERBEAN (SPHERical BEaring ANalysis) has been created 
to simulate the performance of double row spherical roller 
bearings, figure- 1, under a variety of operating conditions. 
Emphasis has been placed on detailing the effects of roller 
skew, roller end to flange contact and change in clearance 
as functions of speed, mounting fits and temperature. 

The complete range of EHD contact considerations has 
been treated in the computation of raceway and flange contact 
detail. A flexible outer ring option is made available for 
simulation of planet bearing performance, where the bearing's 
outer ring is integral with the planet gear. Program capa- 
bilities also include both, steady state and time transient 
temperature mapping of a spherical bearing system. 

Volume I of this three volume report [1]^ describes the models and 
associated mathematics used within SPHERBEAN. The user is 
referred to the material contained therein for formulation 
assumptions and algorithm detail. Volume III [2] describes the 
results of a series of full scale hardware tests, and demonstrates 
the degree of correlation between performance predicted by SPHERBEAN 


1 Numbers 


in brackets 


refer to references at end of report. 
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and measured data. 

The material present in this, Volume II 
guide the user in the correct and practical 
SPHERBEAN. Input and output, guidelines for 
sample executions are detailed in the matter 
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is structured to 
mplementat ion of 
program use and 
which follows. 
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2 . 0 INPUT DATA 

SPHERBEAN requires the preparation of input data which in 
general describes the bearing geometry, properties of the 
materials used, and specifies the imposed operating conditions. 

SPHERBEAN' S input data preprocessor has been structured 
to impose minimal initial demands on the user. Data is segre- 
gated into sets, or CATEGORIES, which individually address 
specific features of the system addressed. Any data set item 
required by the program falls into one of nine distinctly 
identified categories (Table 1). Category "ROLLR" for example, 
contains all roller geometry data. Category "CAGE" details 
cage information. Items detailing planetary gear parameters 
are entered into input category "GEAR". 

All data required by the basic program are accepted in 
free NAMELIST format and default values are hardcoded to mini- 
mize the demands on user judgement. Special input data, required 
when the program options are used, is in 80 column card image 
format . 

In each CATEGORY, the user has the freedom to specify 
all, part, or none of the data. If an item of data is omitted, 
a default value is assumed. A list of these default values 
is shown in Table 2. Failure to include certain basic data, e.g., 
groove radius, load, speed, etc. results in a diagnostic abort 
message. 
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Data comprising a category are specified in free format, 
with the restrictions that: (1) column one of any card is not 
used and (2) all items of data in any CATEGORY are separated 
by commas'. No specific sequence of data is required within 
each CATEGORY, A minimum of two cards is needed to specify a 
complete set of data within a CATEGORY. 

As an example, consider CATEGORY nine, "LIFE. " Here, 
surface finish and fatigue life data are described. Five of 


the required items 

are : 

1 . 

RMSROL - 

The RMS surface roughness of the 
rollers . 

2. 

RMSIR - 

The RMS surface roughness of the 
inner raceway. 

3. 

RMSOR - 

The RMS surface roughness of the 
outer raceway. 

4. 

CIR 

User supplied life multiplier for 
the inner raceway. 

S. 

COR 

User supplied life multiplier for 
the outer raceway. 

An 

example of 

free format input data is illustrated for 


CATEGORY LIFE in Figure 2. The first card contains a dollar 


^ On some computers the commas must follow the last signifi- 
cant digit of an input variable. It is suggested that this 
restriction be observed to avoid the inconvenience caused 
by compiler peculiarities. 
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sign ($) 3 in column 2 followed immediately by the CATEGORY name 
LIFE. The second and third card are used to specify values 
of input data. Note that free format is used throughout and 
that all data items are separated by a comma. The third card 
contains a Dollar Sign ($) in column 2 and the word END in 
columns 3 through 5 signifying the end of data for the category 
Note that column 1 of any card was never used to specify data. 

Each data category is. used to describe a particular aspect 
of either program use the bearing configuration. Categories 
in turn, must be arranged in the sequence noted in Table 1 be- 
fore they can be used to transmit data to the program. 

The following paragraphs will list, in their proper order, 

all categories and the data they contain. In certain cases 

where, at the user's option, categories can be omitted from 

the set, these options are made clear. Likewise, it is also 

indicated when a category must be included, regardless of 

execution options. If the user wishes to omit a category of 

data, he must still include the two cards: 

$ CATEGORY NAME 
$END 


5 Different computers may allow or require a different symbol. 
For example, an IBM 370/158 uses an ampersand (S); a UN I VAC 
1108 uses the dollar sign ($). 
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CATEGORY 1 - PROGRAM LOGIC 
CATEGORY NAME: LOGIC 

CATEGORY DESCRIPTION 

Within this category, the user is permitted to 
specify values for logic used in a given program 
execution. User provided values dictate the program 
options. 

All variables in this category are "logical,” 
and have either of two values, .TRUE, or .FALSE, e.g., 
SYMY = .TRUE.. 

In many cases, additional data will be required 
from the user as a consequence of selecting a speci- 
fic program option. Descriptions of extra input 
data may be found in the Section "SPECIAL INPUT DATA," 
starting on page 28. 


DATA ITEMS: 


DEFAULT : 

FLGIR 

- Used to identify a spherical roller 
bearing having a flanged inner ring. 

F 

PLTROL 

- Allows the user to obtain a line 
printer plot of the roller profile 
along its axial effective length. 

F 

ECHO 

- Allows the user to echo check the 
input data. This option involves 
routines which print the data imme- 
diately after it has been read. 

T 
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DATA ITEMS: DEFAULT 

SPLTCG - Used to execute the split, two piece F 

cage analysis portion of the program. 

The following restrictions must be ob- 
served when employing the split cage 
analysis: 

• use only for pure axial or com- 
bined loading. Under pure radial 
loading, each row of the two row 
bearing is assumed to behave iden- 
tically. 

• use only at solution level 3 or 4 
(LEVEL = 3 or LEVEL = 4 in Cate- 
gory ,, SOLVE M ) 

• do not use with PLANET = .TRUE., 
since each row of a two row plane- 
tary bearing is assumed to behave 
identically. 

PLANET - Used to indicate that the bearing to F 

be simulated is employed in a planetary 
gear application. If .TRUE., the user 
must supply the data items for category 
5, ’’GEAR ” (’’GEAR” described on page 18 ). 

METRIC - If .TRUE. , input data items are being T 

supplied in Metric (M - K - S) units; if 
.FALSE., English (F - P - S) units are 
used. Program output is dimensionally 
consistent with user input. 

Special Program Option Logic : 

The following five data items allow the user to in- 

7 
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yoke certain program options. All items require the 
user to include additional data if a value other than 
the default value is specified (see "SPECIAL INPUT 
DATA"). 

DATA ITEMS: DEFAULT 

FITS - If .TRUE., the bearing's operating diam- F 
etral clearance is computed as a function 
of temperature, speed and initial mount- 
ing fit pressure. If .FALSE., the clear- 
ance change analysis is not executed and 
the user input diametral clearance ("DCL" 
in Category 4) is used in calculating 
bearing details. 

MPROP - Allows the user to input material prop- F 
erties for the rings, rollers, shaft 
and housing. 

SYMY - A .TRUE, value indicates that the roller T 
is symmetric about its y-axis. If SYMY 
is .FALSE., the non-symmetric roller pro- 
file must be read in. 

EVSLIC - The program uses a slicing technique to T 
compute the roller-raceway loads. By 
default all slices are of equal width. 

In many cases it may be advantageous to 
specify slices of unequal width. Un- 
equal slice widths specifying EVSLIC = 

.FALSE. . 

THERM - A .TRUE, value allows the user to use F 

either the steady state or time trans- 
ient temperature calculating routines. 

8 
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CATEGORY 2 - SOLUTION CONTROL PARAMETERS 
CATEGORY NAME: SOLVE 

CATEGORY DESCRIPTION 


Computer program SPHERBEAN uses a modified Newton- 
Raphson iterative scheme to compute values for the 
governing equilibrium equation set. The user may wish 
to override existing solution control parameters. 

Those which are permitted as input by the user are 
contained within this category. 


DATA ITEMS: 
TOL 


NITS 


DEFAULT : 

Convergence criterion used to halt 05 

the iteration procedure. Solution 

is said to be obtained when 

N 7 1/2 

[(Z EQf)/N] < TOL 
i=l 1 

where EQ^ = residue for ith equation 

N = number of equations to be solved 

Maximum number of iterations to be 20 

used in the Newton-Raphson iteration 
scheme. 


LEVEL - Solution control level: 


1 


LEVEL = 1: Inner ring and roller equil- 
ibrium is mutually satisfied 
through consideration of the elas- 
tic (Hertzian) contact loads at 
all roller-race contacts. Roll- 
ers are permitted to translate 
radially, and rotate about their 
z -axis (roller tilt). The inner 
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DATA ITEMS: 

ring may translate along each 
of its three coordinate axes. 
This execution level will per- 
mit economic calculation of 
bearing life. Lubricant and 
friction related loads are not 
considered; epicyclic roller 
speeds are assumed. 

LEVEL = 2: Mutual equilibrium of rollers 

and inner ring is sought as 
in Level 1, then used to eval- 
uate the lubricant related 
forces (EHD, HD, drag, and 
cage pocket). This level of 
execution permits an economic 
estimation of bearing gener- 
ated heat. Epicyclic roller 
speeds are assumed. 

LEVEL = 3: Two solution algorithms are 

used to evaluate the bearing 
performance. The first seeks 
mutual equilibrium of inner 
ring and rollers by varying 
the inner ring position and 
roller radial, axial and tilt 
positions. The second algorithm 
seeks mutual equilibrium of the 
rollers and cage by varying the 
cage position and roller skew, 
tilt and rotational speeds. 
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DATA ITEMS: DEFAULT 

This level of execution 
provides an economic method 
of evaluating bearing perfor- 
mance when there is little 
or no rolling element skid- 
ding and can be used to pre- 
dict the existence of such 
skidding, but not its abso- 
lute magnitude. Note that 
at this level of execution, 
the inner ring reactive load 
may not exactly equilibrate 
with the applied loading. 

LEVEL = 4: The two algorithms used 

in Level 3 are repeated until 
mutually satisfied. This 
assures equilibrium of the 
inner ring with applied 
loading. This execution level 
can be used to predict the 
presence and magnitude of 
rolling element skidding. 

PRINT - Debug output print flag. Allows the 

user to see calculated results , in English 
units, at intermediate steps of solution. 

The PRINT array may be input with the 
following logical variables: 

PRINT (1)=.TRUE. - Displays algorithm number, F 

number of equations to be 
solved, and variables and 

11 
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DATA ITEMS: DEFAULT: 

equations included in 
partition of field equation 
set. 

PRINT(2)=.TRUE. - Elastohydrodynamic (EHD) fric- F 

tion coefficients are printed. 

PRINT (3)=. TRUE. - Roller centrifugal forces are F 

printed. 

PRINT (4) =. TRUE. - Roller gyro-moments are printed. F 


PRINT ( 5) = . TRUE. - Print the corrections of the F 

variable values as calculated 
in the NR equation solver. 

PRINT ( 6) = . TRUE . - Print the matrix of partial F 

derivatives (caution: large 
volume output can be expected) 

PRINT ( 7) =. TRUE. - Print the X and EQ values F 

about which partial deriva- 
tives are computed. 

PRINT(8)=.TRUE. - Print the current NR loop F 

statistics . 

PRINT(9) = .TRUE. - Displays divergence messages, F 


intermediate equation residues 
(EQ) , roller-raceway loads and 
moments, roller independent (X) 
variables, and the solution 
algorithm number. 

PRINT (10)=. TRUE. - Print cage loads and moments. F 
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DATA ITEMS: DEFAULT 

CAGDOF - Number of degrees of freedom to be 1 

used in the cage analysis. Allowable 
integer values are either 1 or 3. The 
one degree-of -freedom cage permits 
cage rotational movement about its 
longitudinal center line; three 
degrees-of-freedom permit translation 
in the radial plane and rotation 
about the centerline. Note that in 
cases of pure axial loading, the cage 
will behave symmetrically and only 
one degree-of-freedom will be employed. 
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CATEGORY 3 - ROLLER GEOMETRY DATA 
CATEGORY NAME: ROLLR 

CATEGORY DESCRIPTION : 

Within this category the user must describe 
the geometry of the rolling elements within the 
bearing complement. This category is always in- 
cluded. Dimensional units are shown in parentheses, 
i.e. (inches [millimeters] ) , where the term in brackets rep- 
resents "metric" system input; roller geometry is illustrated in Pig. 3. 


DATA ITEMS: 


DEFAULT: 

NS 

- Number of roller raceway slices used 
in the analysis. Note: If the user 
specifies symmetry about the roller 
y-axis, i.e. SYMY=.TRUE. in category 
"LOGIC," then NS is the number of 
slices per symmetric half. If the 
user specifies no symmetry then NS 
is the total number of slices. 

5 

HRD 

- Roller maximum diameter (in[mm]) 

None 

ELO 

- Effective length of the roller to 
outer ring contact 4 (in[mmj) 

None 

ELI 

- Effective length of the roller to 
inner ring contact 4 (in[mmj] 

None 


4 The effective contact length refers to the longest possible length 
which can be used to transmit load between the roller and raceway. 
Typically, this is the roller total length less comer radii. If, 
however, the raceway undercuts are exceptionally large so that the 
track width is less than the roller effective length then the track 
width should be input. 

14 
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DATA ITEMS: 

FLENTH - Roller flat length (in [mm]) 

RCC - Roller crown radius (in [mm]) 

KLUCRN - Roller geometry flag having the follow- 
ing possible values: 

KLUCRN * 1; The roller active profile is either 
fully flat or crowned with a flat. 
Symmetry about the y-axis is assumed 

KLUCRN = 2; The roller active profile is spheri- 
cal with y-symmetry (KLUCRN=2 is 
applicable to the roller shown in 
Figure 3). 


KLUCRN = 3; The roller active profile is symmet- 
ric about the y-axis and will be 
read in by the user. (See "SPECIAL 
INPUT DATA"). 

KLUCRN = 4; The roller active profile is non- 

symmetric about the y-axis and will 
be read in by the user (See "SPECIAL 
INPUT DATA") . 


PHI 1 - Angular location of the first rolling 

element in the first row of a double 
row bearing. In the nomenclature of 
Figure 4, the angle is measured CCW 
positive from the bearing y-axis. In- 
put is in degrees. 

PHI 2 - Angular location of the first rolling 

element in the second row. 
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None 

None 
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DATA ITEMS: 

NUMROL - Number of rollers per row. 

REND - Roller end sphere radius Iin[mm]) 

RSREX1 - The x-coordinate of the origin of 
the roller end sphere radius, for 
rollers in row 1 (in[mm]) 

RSREX2 - The x-coordinate of the origin of 
the roller end sphere radius, for 
rollers in row 2 (in[mm]) 

ENDPLY - Roller- flange end play (in(mm]) 
See Figure 5. 
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DEFAULT: 

None 

0 

0 

0 
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CATEGORY 4 - RING DATA AND BEARING APPLIED LOAD 
CATEGORY NAME: RNGEO 

CATEGORY DESCRIPTION : 

Within this category the user must specify 
inner and outer ring geometry, inner ring speed 
and operating load This data set must always be in- 
cluded. (See Figures S and o). 


DATA ITEMS: DEFAULT: 

GROR - Outer ring groove radius (in[mm]). 0 

GRIR - Inner ring groove radius (in [mm]). 0 

DCL - Bearing diametral clearance (infminj). (t 

ALPHA - Absolute value of contact angle (radians) 0 

SPDI R - Inner ring rotational speed (rpm)^ 0 

FLGANG - Flange angle (radians) 0 

5 

1 ■' X - Axial load in X direction (lbs[N|) 0 

S 

FY - Radial load in Y direction (lbs[Nj) 0 

r 

FZ - Radial load in Z direction (lbs[N])° 0 


These data items can be omitted if simulating a planetary bear- 
ing, i.e. PLANET - .TRUE, in Category 1, "Program Logic." Lead 
and speed information are entered in Category 5, "Planet Bearing 
Data." 


!Li i 11 ]'. i V. \::A bLiWlCh; 




17 


AT81DU07 


CATEGORY 5 - PLANETARY BEARING DATA 
CATEGORY NAME: GEAR 

CATEGORY DESCRIPTION : 

SPHERBEAN has the capability to analyze double 
row, spherical planet bearings (Figure 7). In this 
application, the bearing outer ring is integral with a 
gear, and load is transmitted to the bearing through 
a pair of diametrally opposed gear meshes. Figure 8. 
SPHERBEAN will also compute the effects of carrier 
motion on bearing performance by consideration of 
centrifugal forces generated by the kinematics of 
planetary motion, Figure 9. This data set must be 
included when PLANET * .TRUE. 

DATA ITEMS: DEFAULT: 


NMESH - Number of gear tooth pairs in simul- 1 

taneous contact at each mesh. Values 
of 1, 2 or 3 are permitted. 

TANGLE - Angle between teeth on planet gear; 0 

TANGLE * 2n/Number of teeth on pianct 
gear (radians). 

SPDCAR - Carrier rotational speed (rpm). (J 

SPPR1 - Outer ring rotational speed. If carrier 0 
speed is input as non-zero, then SPDR1 


must be input with a negative sign to main- 
tain proper sign conventions (rpm). 

NPL - Number of planet bearings in stage 1 
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DATA ITEMS: DEFAULT : 

RCAR - Radius from sun gear center to 0 

post, Figure 9 (in [mm]). 

XI - Cross sectional moment of intertia 0 

of planet bearing outer ring, Figure 
10. Default is rigid ring (in 4 [mm 4 ]). 

RNUET - . Radius, from bearing centerline, to 0 

neutral axis of the planet bearing 
outer ring, Figure 10 (in[mm]). 

R1 , R2 , R3 - Radial component of gear tooth con- 
tact load at each of up to 3 teeth 
in contact, Figures 11 and 12. 

Number of non-zero input values for 
load components must equal NMESH. 

Padial loads are input with a neg- 
ative sign to maintain sign conven- 
tion, Volume I [1] of this report 
(lbs [N] ) . 

Tl, T2, T3 - Tangential component of gear tooth 0 
contact load at each of up to 3 
teeth in contact, Figure 11 (lbs(N]). 

XM1, XM2, XM3 - Moment component of gear tooth con- 0 
tact load at up to 3 teeth in contact, 

Figure 11 (in-lbs [NmmJ ) . Note that 
XMl-Tl(Rp-RNUET) . 

DENS - Density of the planet gear bearing .28 [.007750] 

outer ring (lbs/in** [KG/cm 3 ] ) . 

XMASS - Weight of the planet gear bearing 0 

outer ring (lbs[N]). 

£ - Modulus of Elasticity 29,500,000 

(lbs/in 2 tN/mm 2 ]) . [204,083] 


SKI HOHNOK JGYSLkVICLS 

SKI INtXJSlMtS.IISK: 


19 



- 41 fli>C - 1 - R100 


ST81D007 


CATEGORY 6 - CAGE DESCRIPTION 
CATEGORY NAME: CAGE 

CATEGORY DESCRIPTION : 

The data items contained within this category are 
used to describe cage geometry. This set of data must 
always be included. (See Figure 13). 

DATA ITEMS: DEFAULT: 

ITYPE - Cage type flag 0 

ITYPE = 0 ; the cage is rolling element riding 

ITYPE = 1 ; the cage is outer ring- land riding 

ITYPE = 2 ; the cage is inner ring- land riding 


RADCLR - Pocket radial clearance (in [mm]) 0 

AXCLR - Pocket axial clearance (in[mm]) 0 

WEBT - Web thickness (in [mm]) 0 

NRAILS - Total number of cage guide rails (in- 1 

eluded for reference only, not used in 
program) . 

A - Inside radius of split interface for 0 

a two piece cage, Figure 14 (in [mm]). 

B - Outside radius of split interface for a 0 

two piece cage. Figure 14 (in[mmj). 

XMU - Coefficient of friction between split .07 

halves of a split cage. 


The following data is included if the cage is land 
riding (ITYPE = 1 or ITYPE = 2): 
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DATA ITEMS: 


RAILD - Equivalent rail land diameter : 

RAILD = 2a/ (a 3 /b)‘ 375 


Where 


and 


NRAILS 


a = u 
i = 1 


i = NRAILS 


b - l 
i = 1 




(r-i) A 


r = the rail land radius 

of the i-th rail (in [mm]). 

£ = the width of the i-th 
rail (in [mm] ) . 

RAILW - Equivalent single rail width^: 

RAILW = (a 3 /b)* 125 

where a and b are defined as above, in 
definition of variable RAILD. 


DEFAULT 

0 


RCLR - Rail land diametral clearance (in [mm]). 


If performing a split cage analysis, i.e. SPLTCG * .TRUE, 
in category LOGIC, these data items are computed for only 
one of the split cage pieces. 
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CATEGORY 7 - LUBRICATION DATA 
CATEGORY NAME: LUBE 

CATEGORY DESCRIPTION : 

Within this category the user specifies lubri- 
cant properties and other data which relate directly 
to the lubricant or to the definition of friction 
related processes. This category should always be 
included. . 


DATA ITEMS: 


DEFAULT: 


NCODE - The user is given the option to 

specify particular lubricant prop- 
ties or simply select a value of 1 
through 4 for NCODE. Specific val- 
ues of NCODE and associated lubri- 
cant properties are shown in Table 
3. The user may input lubricant 
properties not in Table 3 by spec- 
ifying NCODE = 0, and including the 
five specific lubricant properties 
on page 61. 

ZTO, ZTI, ZTFL- Lubricant replenishment layer 

J n | 

thickness' at the outer raceway, 1 

inner raceway and the flange, 


20. X 10' u . 

[5.08 X 10' 4 ] 


At the present time the magnitudes of the inner and outer replenishment layer 
thicknesses have not been correlated with flow rate, particular lubricants or 
bearing speed. The user is required to establish proper values of the replen- 
ishment layer thickness. The following guidelines are suggested: 

•To avoid starvation ZTO, ZTI, ZTFL should be 1 or 2 times the 
EHD film thickness 

•Centrifugal loading suggests ZTO be larger than ZTI. 
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DATA ITEMS: DEFAULT: 

respectively (in/mm). 

Note: if ZTO = ZTI = ZTFL = 

0.0 dry friction is assumed. 


FRK - Lubricant pseudo friction co- .07 

efficient, used in the Allen [ 3 ] 
traction model. Typical values 
lie in the range 0.5 ^ FRK ^ .08. 

AKN - SPHERBEAN uses a model developed 18.2 

by Loewenthal [ 4 ] to compute EHD 
film thickness. The term AKN, 
the lubricant film thickness co- 
efficient, appears in that equa- 
tion. Typical values are 
18.2 $ AKN <50. 


XMURC - Dry coefficient of friction at .075 

race contacts. If ZTO-ZTI=0, 

XMURC is applied. 

XMUFL - Dry coefficient of friction at the .075 

flange contact. If ZTFL = 0, 

XMUFL is applied. 

XMUCG - Coulomb friction coefficient .075 

used at the cage pocket-rolling 
element contact. If ZTO=ZTI=0, 

XMUCG is applied. 


XCAV - Percent of lubricant occupying 1, 

the bearing cavity®. O.s? XCAV< 100. 

8 7 7 

As with replenishment layer thickness, the amount of free lubricant should 
be correlated with the operating parameters. At this time, such correla- 
tions do not exist. XCAV values of less than 5 percent are recommended. 
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DATA ITEMS: 


DEFAULT : 


BRKPT - Limiting value of slide-to- 0.005 

roll ratio used in calculating 
friction forces with the Allen [3] 
traction model when the critical 
shear stress is exceeded. 

The following data items must be included if NCODE was specif i 
as zero. 


DATA ITEMS: 


DEFAULT: 

RH060 

3 

Density of lubricant (gm/cm ) at 
60°F. 

None 

G 

Thermal coefficient of expansion 
(1°F) . 

None 

VIS100 

Viscosity of lubricant (centi- 
stokes) at 100°F. 

None 

VIS210 

Viscosity of lubricant (centi- 
stokes) at 210°F. 

None 

COND 

Thermal conductivity (BTU/hr- 
ft-°F) 

None 
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CATEGORY 8 - OPERATING TEMPERATURES 
CATEGORY NAME: TEMPS 

CATEGORY DESCRIPTION : 

Operating temperatures are defined in this 
category. Outer race, inner race and flange temp- 
eratures are used to evaluate the properties of the 
specified lubricant at these locations. Bulk temp- 
erature (BULK) is used to evaluate the properties of 
the lubricant contained in the free space of the bearing 
cavity. This information is subsequently used in the 
calculation of the viscous drag force acting upon the 
rolling elements. Data for this category must be 
included unless performing a thermal analysis, i.e. 

THERM = .TRUE, in category LOGIC. Temperatures are 
specified in either °F or °C. 

DATA ITEMS: (See Figure 15) DEFAULT: 

RING 1 1 - Outer ring temperature for row no. 1. 212. [100.] 

RING12 - Outer ring temperature for row no. 2. 212. [100.] 

RING21 - Inner ring temperature for row no. 1. 212. [100.] 

RING22 - Inner ring temperature for row no. 2. 212. [100.] 

BULK - Average temperature of lubricant in 212. [100.] 
bearing cavity. 

FLG1 - Inner ring flange temperature for 212. [100.] 

row no . 1 . 
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DATA ITEMS: 

FLG2 - Inner ring flange temperature for 
row no. 2. 

SHAFT - Shaft temperature 

HOUSE - Housing temperature 

R0LL1 - Temperature of the rollers for row 
no. 1 . 

ROLL 2 - Temperature of the rollers for row 

no . 2 . 

RACE1 1 - Outer race temperature for row no. 1. 

RACE12 - Outer race temperature for row no. 2. 

RACE 21 - Inner race temperature for row no. 1. 

RACE 2 2 - Inner race temperature for row no. 2 , 

CAGE1 - Cage temperature for row no. 1. 

CAGE2 - Cage temperature for row no. 2. 


DEFAULT: 
212. [100.. 1 

212 . 1100 .] 

212 . 1100 .] 

212 . 1100 .] 

212 . [ 100 .] 

212 . [ 100 . ] 
212 . [ 100 . ] 
212 . [ 100 . ] 
212 . [ 100 . 1 
212 . [ 100 . ] 
212. [100 . ] 
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CATEGORY 9 - SURFACE FINISH AND FATIGUE LIFE DATA 
CATEGORY NAME: LIFE 

CATEGORY DESCRIPTION : 

Self-explanatory. This category must always 
be included. 

DATA ITEMS: DEFAULT: 


RMSROL 

- The RMS surface roughness 
of the roller (in [microns] ) 

8x10 _6 [.20] 

RMSIR 

- The RMS surface roughness of 
the inner ring (in [microns] ) 

IOxIO 6 [ * 25] 

RMSOR 

- The RMS surface roughness of 
of the outer ring. (in [microns ] ) 

10xl0' 6 [. 25] 

RMSFL 

- The RMS surface roughness of 
the flange (in[microns] ) 

10xl0' 6 [. 25] 

RMS RE 

- The RMS surface roughness of 
the roller end (in [microns ] ) 

10xl0' 6 [.25] 

CIR 

q 

- Life correction factor for the 
inner ring. 

1 . 

COR 

9 

- Life correction factor for the 
outer ring. 

1 . 


3 

The numbers input for CIR and COR are used to account for improved mater- 
ials by multiplying the raceway fatigue lives as calculated by Lundberg- 
Palmgren methods. Typical life factor values' for modem steels are in the 
range of 2 to 3. 

In the ASME Publication Life Adjustment Factors for Ball and Roller Bearings , 
the Material Factor D and the Material Process Factor E should be used 
multiplicatively as inputs for CIR and COR. 
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3.0 SPECIAL INPUT DATA 

The user activates the extended program capabilities by 
invoking a maximum of up to six program options. Logic used to 
activate a given option was presented in the previous section, 
and is summarized in Table 4. 

The user may employ as many options as necessary in a 
single program execution. The only restrictions are found in 
the input data sequence shown in Table 4 and the use of the 
specific 80 column card formats detailed in Appendix A. 

All options require the user to specify additional 
information. This data follows immediately after the basic 
categorized data. 

! 

j 

i 

I 


u. 
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3.1 OPTION 1: COMPUTATION OF BEARING OPERATING 
DIAMETRAL CLEARANCE 

This data should be included only if the change in 
bearing diametral clearance, from "off the shelf" to operating 
value is to be calculated. This option is invoked by speci- 
fying FITS » .TRUE, in category LOGIC. The user must supply 
effective widths and radii of the bearing components, shaft 
and housing. Positive values should be used for interference 
fits and negative for clearance (measured on the radius). All 
data must be given in the cold state (68°F). 

Card formats for data input are shown in Figure A1 . Input 
data items are defined in Figure A1A. 
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3.2 OPTION 2; USER SPECIFIED MATERIAL PROPERTIES 

The user may specify the material properties of bearing 
components by setting the logical variable MPROP = .TRUE.. 
Appropriate properties arc input according to the card format 
shown in Figure A2. Unspecified variables (i.e. those set 
to ’,ero or left blank) are assigned the following values: 


Modulus of Elasticity 

Poisson's Ratio 

Coefficient of Thermal 
Expansion 

Density 


204085 N/mm 2 
( 2 . 9 59988xl0 ? psi ) 

0.3 

1.124 x 10" 5 per °C 
(6.744 x 10’ () per °F) 

7.806 gm/cm 3 
(0.28 lb/in 3 ) 
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3.3 OPTION 3: USER SPECIFIED SLICE WIDTHS 

SPHERBEAN uses the slicing technique (see Volume 1 [1] of 
this report) to compute the roller to raceway contact traction 
forces. Ordinarily, the slice widths are assumed equal, how- 
ever, in some instances it may be advantageous to specify 
their individual and varying extent. 

This option may be invoked by specifying SYMY = .TRUE, 
and EVSLIC = .FALSE.. Since symmetry is assumed, the slice 
widths need only be input for a symmetric half of the roller. 
The numbering scheme used is one where the first slice is en- 
countered at the roller centerline, the last at the roller 
end, Figure 16. The card format shown in Figure A3 is used 
for input data. 
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3 . 4 OPTION 4: USER INPUT SYMMETRIC ROLLER GEOMETRY 

The user may input values for roller radii at specific 
locations along the roller effective length. The input 
sequence is shown in Figure 16. 

Note that the number of radii to be read in is equal to 
the number of slices (per symmetric half) plus one. Card 
formats are given in Figure A4. Logic used is SYMY = .TRUK, 
and KLUCRN = 3. 


i. 

r 
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3.5 OPTION 5: USER SPECIFIED, COMPLETELY VARIABLE, 

ROLLER GEOMETRY 

The user can specify the complete detail of roller geom- 
etry by using the options: SYMY * .FALSE, and KLUCRN * 4. 

When using this option the user must specify roller 
radii and slice widths across both the outer and inner ring 
effective lengths. Note that the number of radii input cor- 
responds to the total number of slices plus one, Figure 17. 

Card formats for data input are shown in Figure A5. 
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3.6 OPTION 6: TEMPERATURE CALCULATIONS 

SPHERBEAN may be used to compute either the time transient 
or steady state temperature distribution within a system de- 
fined by the bearing and its environment. Logic used requires 
THERM = .TRUE.. 

The temperature portion of SPHERBEAN is designed to pro- 
duce temperature maps for an axisymmetric mechanical system of 
any geometrical shape. The physical system is first approxi- 
mated by an equivalent "nodal network" which consists of a 
number of elements having simple geometric shapes. Each 
element is then represented by a node point characterized by 
a mass, surface area, and having either a known or an unknown 
temperature. The environment surrounding the mechanical system 
is also represented by one or more nodes. With the node points 
selected, heat balance equations are formulated by the program 
for the nodes of unknown temperature. These equations become 
non-linear when there is radiation between two or more of the 
node points considered. 

The success of the approach depends largely on the realist- 
ic physical subdivision of the system. If the subdivision is 
too fine, there will be a large number of equations to be 
solved, and execution cost will rise. If the subdivision is 
too crude, the results are likely to be inaccurate. 

34 

SKF TECHNOLOGY SERVICES 

Ski- INDUS R^ltb. iNc 


F-4185C-1-R100 


AT81D007 


The present thermal simulation is restricted to the treat- 
ment of axially symmetric physical systems. Bearing rings for 
example, fall into this category and can be represented by an 
element of uniform temperature. For a component or module 
which is not axially symmetric, the user must represent it with 
an equivalent axially symmetric element of approximately the 
same surface area and material volume. 

With input data prepared as described in the following 
paragraphs, SPHERBEAN will solve the heat-balance equations 
for either the steady state or the time transient conditions 
and produce temperature maps for the physical system. 

DESCRIPTION OF INPUT DATA FOR TEMPERATURE CALCULATIONS 

Card formats for data input are listed in Figures A6 . 

Card 1 

Card 1 is a control card and contains input for both 
steady state and transient thermal analyses. It is not in- 
tended however, that both analyses be executed with the same 
run. 

Item 1: Highest Node Number (M) . The temperature nodes 

must be numbered consecutively from one (1) to the highest 
node number. The highest node number must not exceed one 
hundred (100). 
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Item 2: Number of Unknown Temperature Nodes (N) . It is 

required that all nodes with unknown temperatures be assigned 
the lowest node numbers. The nodes which have known tempera- 
tures are assigned highest numbers. 

Ite m 3 : Common Initial Temperature (TEMI)°C: The 

temperature solution iteration scheme requires a starting point, 
i.e,, guesses of the equilibrium temperatures, however, when 
a node is not given a specific initial temperature, the temp- 
erature specified as Item 3 of Card 1 is assigned. 

Item 4 : Punch Flag (IPUNCH): If the Punch Flag is not 

zero (0) or blank, the system steady state equilibrium temp- 
eratures, along with the respective node numbers, will be 
punched according to the format of Card 2. This option is 
useful if, for instance, the user makes a steady state run 
with lubrication, and then wishes to use the resultant temper- 
ature as the initiation point for a transient dry friction run 
in order to assess the thermal consequence of lubricant flow 
termination. 

Item 5 : "Output Flag" (IUB). If the "Output Flag" is 

not zero the bearing program output and a temperature map will 
always be printed after each call to the bearing solution 
scheme. 
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Item 6 : ''Maximum Number of Calls to the Bearing Program" 

(IT1). IT1 is the limit on the number of Thermal -Bearing 
iterations, i.e., the external temperature equilibrium calcu- 
lation^ The user must input a non-zero integer such as 5 or 
10 for SPHERBEAN to iterate to an equilibrium condition. 

If IT1 is left blank or set to zero (0), the number of 
iterations is set to three (3). 

Ill also serves as a limit on the transient temperature 
solution scheme, by limiting the number of times the bearing 
solution scheme is called. Each call to the bearing scheme 
will input a new set of bearing heats to the transient temper- 
ature scheme until a steady state condition is approached or 
until the transient solution time-up limit is reached. 

Item 7 ; "Absolute Accuracy of Temperatures for the 
External Thermal Solution” (EP1) . In the steady state thermal 
solution scheme, each calculation of system temperatures 
occurs after a call to the bearing scheme which produces bearin 
generated heats. After the system temperatures have been 
calculated for each iteration, using the internal temperature 
solution scheme, each node temperature is checked against the 

See Figure 14 (Volume I [1] of this report) for details 
of the thermal solution. 
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nodal temperature at the previous iteration. 

If i^EPl for all nodes i then equilibrium 

has been achieved and the iteration process stops. 

Item 8 : "Iteration Limit for the Internal Thermal 

Solution" (IT2). After each call to the bearing analysis, the 
program will use a Newton-Raphson iteration method to solve 
the heat balance equations for the steady state equilibrium 
temperatures, based on the set of bearing heat generation rates. 
If IT2 is left blank or set to zero (0), the number of Newton- 
Raphson iterations is limited to twenty (20). 

\ 

Item 9 : "Accuracy for Internal Thermal Solution" (EP2) . 

The use of EP2, the convergence tolerance for the Newton- 
Raphson procedure described in Item 8, is explained in Volume I. 
If EP2 is left blank or set to zero (0), a default value of 
0.01 is used. 

Item 10 : "Starting Time" (START) is a time at which the 

transient solution begins, T ; usually set to zero (0). 

s 

Item 11 ; "Stopping Time" (STOP) is the time in seconds 
at which the transient solution terminates, T £ . The transient 
solution will generate a history of the system performance which 

^See Figure 14, Volume I [1]. 
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will encompass a total elapsed time of 

(Tf - T g ) seconds. 

Item 12 : "Calculation Time Step" (STEPIN) . The transient 

internal solution scheme solves the system of equations (see 
Volume I) : 



AT- STEPIN 

The user may specify STEPIN. If left blank or set to 
zero (0), SPHERBEAN calculates an appropriate value for 
STEPIN using the procedure described in [5]. 

Item 13 : "Time Interval Between Printed Temperature Maps" 

(TTIME) seconds. The user must specify the length of time which 
will elapse between each printing of the temperature map. The 
interval will always be at least as large as the "calculation 

time step" (STEPIN) . 

Item 14 : "Time Interval Between Calls of the Bearing 

Program" (BTIME) . BTIME will always have a value larger than 
or equal to (STEPIN) even if the user inadvertently inputs a 
shorter interval. Computational time savings result if BTIME 
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is greater than STEPIN, however, accuracy might be lost. 

Card 2 

In the steady state analysis this card is used to input 
initial guesses of individual nodal temperatures for unknown 
nodes as well as the constant temperatures for known nodes, 
such as ambient air and/or an oil sump. 

In the transient analysis. Card 2 is used to input the 
nodal temperatures of all nodes at time = T g , i.e., at the 
initiation of the transient solution. 

Card 3 

With this card, node numbers are assigned to the components 
of the bearing. With this information the proper system temper- 
atures are carried into the respective bearing analysis. The 
inner race and inner ring node numbers may or may not be the 
same at the user’s discretion. Similarly, the outer race and 
outer ring node numbers may or may not be the same. 

Card 4 

The bearing analysis accounts for frictional heat gener- 
ated at different locations in the bearing, i.e., the inner 
race, the outer race, between the cage rail and ring land, the 
bulk lubricant due to drag and at the flange. The heat 
generated at the cage-rolling element contact is added to the 
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bulk lubricant. This card allows the heat generated to be 
distributed equally to two nodes. For instance, the heat 
generated at the inner race-rolling element contact should be 
distributed half to the rolling element and half to the inner 
race. The heat developed between the cage and inner ring land 
may be distributed half to the inner ring and half to the 
cage if a cage node has been defined, otherwise, half to the 
bulk lubricant. 

Card 5 

This card specifies the node numbers and the heat generation 
rate at those nodes. This card is used to specify where heat 
is generated at a constant rate such as at rubbing seals or 
gear contacts. 

Card 6 

This card type is used to input the numerical values of 
the various heat transfer coefficients which appear in the 
equations for heat transfer by conductivity , free convection , 
forced convection , radiation and fluid flow . Up to ten coeffi- 
cients of each type may be used. Separate values of each type 
of coefficient are assigned an index number via card 6. When 
describing heat flow paths (Card 7 below) it is necessary only 
to list the index number by which heat transfers between node 
pairs . 
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Indices 1-10 are reserved for the conduction coefficient 
X, 11-20 for the free convection parameters, 21-30 for forced 
convection , 31-40 for emissivity and 41-50 for fluid flow 
(product of specific heat, density and volume flow rate). 

As an example, for heat transfer by conduction with 
coefficient X of 53.7 watts/M°C one could prepare a card 6 
with the digit 1 punched in column 10 and the value 53.7 
punched in the field corresponding to card columns 11-20. If 
a conduction coefficient of 46.7 were applicable for certain 
other nodes in the system one could punch an additional card 
assigning index No. 2 to the value X = 46.7 by punching a "2" 
in card column 10 and 46.7 anywhere within card columns 11-20. 

Rather than providing constant forced convection coeffi- 
cients, these coefficients can be calculated by the program 
in one of three ways. If this calculation option is 
exercised a pair of cards is used in place of a single card 
containing a fixed value of a. The contents of the pair 
depends upon which of the three methods are used. 

Option 1) a is independent of temperature but 
is caluclated as a function of the 
Nusselt number which in turn is a 
function of the Reynolds number R^, 
the Prandtl number P as follows, 

(cf. 5 ): 
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Option 2) 


where A is the lubricant conductivity, 
oil 

L is a characteristic length (with the 
units of meters) and K, a and b are con- 
stants . 

a is a function only of fluid dynamic 
viscosity and viscosity is temperature 
dependent. 

a = cn^ 


where c and d are constants. 

Option 3) a is again a function of the Nusselt, 
Reynolds and Prandtl numbers as in 
Option 1, but the viscosity is temper- 
ature dependent. 


Appendix B has been included to aid the user in data 


preparation and calculation of heat transfer coefficients. 


Card 7 

This card defines the heat flow paths between pairs of 
nodes. Hvery node must be connected to at least one other 


node, i.e., two or more independent node systems may not be 


solved with a single program execution. 

The calculation of heat transfer areas is based on 


lengths, and input using card 7. Additionally, the type 
of surface for which the area is being calculated is indicated 
by the sign assigned to the heat transfer coefficient index. 
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If the surface is cylindrical or circular, the index should he 
positive, if the surface is rectangular, the index provided should 
be a negative integer. 

In the case of radiation between concentric axially symmetric 
bodies, is the radius of the larger body. For radiation be- 
tween two parallel flat surfaces or for conduction between nodes, 

Lj is the distance between them. 

Fluid flow heat transfer computation accounts for the energy 
which the fluid transports across a node boundary. Surrounding 
a fluid node at which convection is taking place, the temperature 
varies. The nodal temperature which is computed is the average 
of the fluid temperature at the outlet and inlet boundaries. If 
the emerging temperature of the fluid is of interest, it is 
necessary to have a fluid node at the fluid outlet. At this 
auxiliary node, heat transfer occurs by fluid transport only 
and the fluid temperature is considered to remain throughout 
the volume it represents. The true fluid outlet temperature will 
be obtained in this way. 

Conduction of heat through a bearing is controlled by index 
51. The actual heat transfer coefficient which contains a con- 
ductivity, area and a path length term is calculated in the bear- 

1 2 

mg portion of the program. The term is based on a computation 
which averages the outer race and inner race rolling element contacts. 

12 ... 

See Volume I, Appendix C for a detailed description. 
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4.o output DATA 

All data generated by an execution of SPHERBEAN is suitable 
for output on a 132 character width line printer. Output data 
is organized so that user supplied information, such as bearing 
geometry and operating conditions, appears first. Default value 
for unspecified data are also displayed. 

All computed data, representing the results of the 
analysis performed, is presented after the display of user 
supplied data. 

Two examples will be used to illustrate typical program 
output. The first addresses the steady state thermomehcanical 
simulation of a radially loaded 40mm bore spherical roller 
bearing. The second shows typical output obtained when 
simulating planetary gear spherical roller bearing performance. 
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4 . 1 EXAMPLE 1: Simulation of Steady State Bearing 

Performance 

SPHERBEAN was used to simulate steady state thermomechani - 

cal performance of the roller bearing test rig shown in figure 

18. The system incorporates a 40mm bore double row spherical 

roller bearing (Table 5) , a dummy cylindrical roller bearing, hydro- 

1 ^ 

static radial and thrust bearings, shaft, housing and lubrication systems. 

The test rig and bearings are modelled by the 50 nodes 
shown in Figures 19 and 20 using the program's "Option 6". 

Nodal representations of pertinent rig components are identi- 
fied in Table 6. 

The output generated for Example 1 is presented in Appendix C. 
In this appendix, page 1 lists the user defined spherical bearing 
nodes and nodes where bearing heat is generated. Note that each 
source of bearing generated heat splits its output between 
two bearing components. For example, in row #1, heat 
generated at all roller (node 10) to outer race (node 12) con- 
tacts is evenly distributed: 501 to the rollers and 50$ to the 

outer race. 

Constant heat generation sources are also listed on output 
page 1. In this example, nodes 24 and 28 represent lip seals, 

13 

A complete description of the test rig can be found in 
Volume III of this report. [2]. 
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and 6.7 watts were provided at each site to represent the 
sliding power loss. Nodes 32,33 and 44 were used to simulate 
power loss in the test rig slave hearing. The value of 
180 watts (60 watts at each node) was computed using a 
handbook estimate of the slave bearing’s power loss at 3000 rpm. 

Pages 2 through 6 are an organized listing of the user 
specified special input data needed for system temperature 
calculations. This data completely defines the thermal simu- 
lation. 

Page 7 displays the first system temperature map, and 
represents the initial common temperature of 93°C inpuc on 
Card 2. Nodes 4S through 30 have fixed temperatures, and 
represent the system boundary conditions. 

Page 8 lists options in effect during the program exe- 
cution. Roller and ring geometry, as input by the user, is 
also listed. 

Page 9 presents an organized listing of the user defined 
cage geometry, material properties and lubrication data. 

Surface geometry and applied loads are listed on page 10. 

The bearing fatigue life, as well as individual L 1() 
fatigue lives of the outer and inner rings, arc presented on 
page 11. The bearing life represents the statistical combina- 
tion of the two raceway lives. Raceway lives reflect the 
combined effects of: 
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(1) the user input material factors 

(2) lubricant film thickness factors. 

The film thickness to surface roughness ratio is used in the 
calculation of a life reduction factor. Detailed information 
for this calculation is given in [6] and [7]. 

(3) Life modification for material other than 
basic steel [l 3 . 

(4) Life modification for edge loading of 
rollers [1]. 

The lubricant data shown on page 12 is self-explanatory. 

Temperatures at which properties are evaluated correspond to 

the calculated steady state operating conditions. 

Pages 13-14: The roller-raceway contact loads at the 

outer and inner ring contacts are defined in the R coordinate 14 

frame, illustrated in Figure 21. These forces will include 

both elastic and lubricant traction effects for solution 

levels 2, 3, or 4 . Rollers are numbered in ascending order 

beginning with the roller lying on the bearing y-axis and 

proceeding counter clockwise, Figure 4. 

Page 15: The forces applied to the inner ring constitute 

system loading information specified by the user. Radial 

misalignment is included in the program output for future use. 

The calculated "inner ring reactive forces and moments" are 

the resultant of the vector sum of all roller-inner ring contacts. 

T4 — ] 

The R coordinate frame is fully defined in Volume I, Appendix A [Ij. 
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Load values are included with the output to assess the 
degree of convergence at the lower solution 
levels. In this LEVEL = 2 execution, the user specified 
applied load was 3000 lbs in the Y direction. The program 
arrived at a solution of: 

Y = -2999 lbs 
Z = -87.98 lbs 

Note that the theoretically correct solution requires Z = 0 
and the 87.98 lb Z directed residual represents 3 % deviation 
from the applied load. If the user requires a "tighter” solu- 
tion, SPHERBEAN must be executed at LEVEL = 3 or LEVEL = 4. 

Page 16 lists the centrifugal iorce vectors and Z directed 
gyro-moments at each roller location in the bearing complement. 
If a planetary bearing is analyzed, these loads will include 
the effects of planetary motion. The coordinate system used 
is that shown in Figure 21. 

Page 17 summarizes bearing heat generation rates. Hydro- 
dynamic (HD) or "rolling" friction refers to the power loss 
created by churning the lubricant in the low pressure zone 
preceeding the Hertzian portion of the EHD roller to raceway 
contact. Bearing conductance is a measure of the resistance 
to heat flow across the bearing and is used in the computation 
of the steady state system temperatures. 

Page 18: The Hertzian contact stresses represent maxi mum 
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values for the line (roller-raceway) contact. The lubricant 
film thicknesses represent minimum values for the line contact. 
EHD friction coefficient is computed by summing the absolute 
values of the EHD and HD loads at all slices, then dividing 
this value by the elastic (Hertzian) contact load. 

Page 19 lists the starvation coefficient and meniscus 
distance at each roller location. In this version of SPHERBEAN, 
the starvation coefficient is not applied to the computed film 
thickness, and is included in the program output for the user’s 
reference only. The meniscus distance is the distance from 
the contact center to the edge of the resevoir preceeding 
the Hertzian contact. This value is used in the computation 
of rolling (HD) friction power loss (see [1] for a more 
detailed description). 

1 5 

Page 20 displays the sliding velocity distribution 
across the inner and outer ring contacts for roller number one. 

In this execution roller skew was not considered^, therefore 
the x component of sliding is zero. Note that at the inner 
raceway, slice numbers 1, 2, 3, 18, 19 and 20 are nut in con- 
tact with the inner race, and consequently all components of 
sliding velocity at these slices are printed as zero. 

"^Sliding velocity refers to the difference in surface velocities 
at the point of contact between two slices. 

■^Skew is only considered at LEVEL = 3 or 4. 
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Page 21: Epicyclic cage speed is printed along with the 

calculated roller speeds for the user's reference. Epicyclic 
speeds are those speeds assumed by the bearing components in the 
absence of gross slip. These are useful in assessing roller skid. 
Roller rotational speed and roller orbital speed are computed in 
one of two ways. At solution levels one or two, the printed speeds 
are epicyclic. At levels three and four, the speeds are computed 
along with roller displacement independent variables. 

Calculated roller skew and tilt is presented to the user in 
the roller reference frame (R) , and refers to the rotation the 
roller experiences relative to its initial position about its 
Y (skew) and Z (tilt) axis. Cage pocket leads presented on 
page 22 represent the loads experienced by the roller due to 
interaction with the cage web. In row 1, a negative sign indi- 
cates that that roller is pushing the cage. In row 2, a positive 
sign indicates the roller is pushing the cage. 

Pocket loads are presented for solution levels 2, 3, or 4. 

At level two, pocket loads are estimated from rolling element 
Z equilibrium force equation residues. Pocket loads are computed 
such that rolling elements and cage are in equilibium, at levels 
three or four. 

Page 23 lists cage rail and split cage data. Cage centri- 
fugal load refers to the centrifugal force generated by carrier 
motion in a planetary bearing application. 
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Final operating temperatures of the bearing, as calculated 
under the temperature calculating program option, are shown for 
all nodes on page 24. These steady state temperatures are in 
degrees Celcius. 
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4 . 2 EXAMPLE 2: Planetary Bearing Analysis 

SPHERBEAN output, when the program is used to simulate per- 
formance of a planetary gear bearing having a flexible outer ring, is 
similar to output obtained when analyzing mounted bearings. 

Three output pages which will be different are included in 
Appendix D. 

Page 1 (which is the equivalent to page 10 in Appendix C) 
lists planet gear geometry, speed and load data. Note that in this 
example the carrier speed is assumed to be zero, and the radial 
component of gear tooth load has been (correctly) input with a neg- 
ative sign. These valves are entered in Category (5) . 

Page 2 (equivalent to page 15 in Appendix C) details the 
applied and reactive inner ring loads. Note that the diam- 
etrally opposed tangent and moment tooth loads will cancel 
each other and produce no net load on the inner ring. The 
inner ring or "Post Load" will be the vector sum of tangent 
gear loads and the outer ring centrifugal load due to carrier 
rotation (which in this case is zero). 

Page 3 is only printed when the planetary bearing analysis 

option is used, and lists the deflection of the outer ring at 

individual roller locations. The sign convention employed is 

1 7 

positive outward. 

17 

This sign convention is consistant with the one used tor in- 
put of planet gear tooth radial load. 
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5.0 PROGRAM LIMITATIONS AND GUIDES TO PROGRAM USE 

SPHERBEAN is a design tool, and as with any tool, successful 
use requires awareness of intended applicability and inherent 
limitations. 

LIMITATIONS 

The user must conform to the following geometric and opera- 
ting restrictions: 

(1) The bearing complement may contain no more than 
thirty rollers per row. However, in order to 
execute the program at levels 3 or 4 under combined 
load for double row bearings having more than 17 
rollers per row, it is necessary to increase the 
size of the C array which stores partial derivatives. 

The array is currently dimensioned to accommodate a 
144 X 144 array of partials [COMMON/LINCO/C(20740) ] .. 

In the event that the array size must be increased, 
the program prints out the error message given 
below indicating the required dimension size and 
execution is stopped. 

NUMBER OF DERIVATIVES HAS EXCEEDED THE 
ALLOWABLE STORAGE OF 20740 

SPHERBEAN NEEDS AT LEAST STORAGE LOCATIONS 

TO DO CALCULATIONS FOR YOU 
DIMENSION LIMITS MUST BE INCREASED* ***** 
EXECUTION STOPS ***** 

The /LINCO/ common appears in subroutines DERIVS, SIMQ, 
and WRITC. It is also necessary to change the value of 
MAXNZ in subroutine DERIVS to correspond to the required 
C array size. 

(2) System thermal models may use no more than 100 nodes. 
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(3) The number of unique heat transfer coefficients is 
limited to: 

• conduction 10 

• free convection 10 

• forced convection 10 

• radiation 10 

• mass transport 10 

(4) The number of heat transfer links is limited to 500. 


GU IDES TO PROGRAM USE 

SPHERBEAN uses a modified Newton-Raphson iterative method 
to compute solution values to the governing field equation set. 

The hard coded (default) values for iteration loop control, con- 
tained in CATEGORY SOLVE, were selected to provide satisfactory 
program performance for a broad class of bearing problems. In 
special instances, where the problem addressed is particularly 
complicated or bearing performance is being checked at an opera- 
ting extreme the default parameters may be inadequate to provide 
a satisfactory solution. In this instance, the user may choose 
to override the default values. The following guidelines are meant 
to assist the user in defining when and what corrective action 
could be taken: 


1. No corrective action need be taken if the 
diagnostic message shown in Table 7-A is 
printed. The message simply indicates that 
a problem was encountered during the solu- 
tion procedure, but was self-corrected. 

2. If the error message in Table 7-A is followed 

by 7-B, this indicates that the desired solution 
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accuracy has not been achieved and results 

1 8 

are not printed. The user can take the 
following corrective action: 

* Simplify the problem - Ex.: If CAGDOF = 3, 
set CAGDOF = 1. 

* Note the value for CURRENT NORM. If it is 
"close" to the value printed for REQUIRED 

1 q 

NORM, the user may input a value of "TOL" 
which equals 1/2 the value printed for the 
CURRENT NORM. An upper limit of TOL = .05 
is suggested. 

* Try a lower solution level. If LEVEL = 3 
or 4, change to LEVEL = 2. 


In some instances the tolerance (NORM) is very close to the 
desired value when the message is printed. If the current 
norm is within 100% of the user specified requirement when 
no further improvement could be made, the program will 
continue the execution. 

See CATEGORY SOLVE. 
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CATEGORY 

NUMBER 

CATEGORY 

NAME 

CATEGORY DESCRIPTION 

INPUT 

DATA 

NECESSARY 

1 

LOGIC 

PROGRAM LOGIC 

NO 

2 

SOLVE 

SOLUTION CONTROL DATA 

NO 

3 

ROLLR 

ROLLER GEOMETRY 

YES 

4 

RNGEO 

RING GEOMETRY/LOADS 

YES 

5 

GEAR 

PLANET GEAR DATA 

NO 

6 

CAGE 

CAGE DESCRIPTION 

YES 

7 

LUBE 

LUBRICANT DATA 

NO 

8 

TEMPS 

TEMPERATURES 

NO 

9 

LIFE 

FATIGUE COMPUTATION 
DATA 

NO 


TABLE 1 

ORGANIZATION OF INPUT DATA CATEGORIES 
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ORGANIZATION OF SPECIAL INPUT DATA 
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Rings: 


t 

i* 


Rollers 


Cage: 


Outer Ring Groove Radius 

1. S898 

in. 

Inner Ring Groove Radius 

1.6923 

in . 

Operating Diametral Clearance 

. C02 

in 

Contact Angle 

14.17 

DEG 

Inner Ring Speed 

3000 

RPM 

Inner Ring Load 

3000 

LBS 

Diameter 

.5122 

in . 

Crown Radius 

1.574 

In. 

End Sphere Radius 

5.552 

in . 

Effective Length 

. 469 

in. 

Number of Rows 

2 


Number of Rollers Per Row 

12 


Pocket Radial Clearance 

.01 

in. 

Pocket Axial Clearance 

0 

in. 

Pocket Web Thickness 

.137 

in. 


(Roller Riding Single Piece Cage) 


i 

TABLE 5 

SPECIFICATIONS OF SPHERICAL ROLLER 
BEARING USED IN EXAMPLE 1 

> o 

I 2 

0: 
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NODE NUMBER 

4 

6 

7 

8 
9 

10 

11 

12 

13 

14 
19 
22 
24 
28 

32 

33 
36 

36 
44 
46 
4 b 


4b 

49 
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DESCRIPTIO N 

SHAjT 

SPHERICAL BEARING - INNER RING ROW #1 
SPHERICAL BEARING - INNER RING ROW if 2 
SPHERICAL BEARING - CAGE ROW #1 
SPHERICAL BEARING - CAGE ROW #2 
SPHERICAL BEARING - ROLLERS ROW #1 
SPHERICAL BEARING - ROLLERS ROW # 2 
SPHERICAL BEARING - OUTER RING ROW #1 
SPHERICAL BEARING - OUTER RING ROW it 2 
HOUSING 

COMPARTMENT AIR 
COMPARTMENT AIR 
LABYRINTH SEAL if 1 
LABYRINTH SEAL if 2 

DUMMY CYLINDRICAL BEARING - ROLLERS 
DUMMY CYLINDRICAL BEARING - INNER RING 

LUBRICANT IN DUMMY CYLINDRICAL BEARING 
CAVITY 

LUBRICANT IN SPHERICAL BEARING CAVITY 

DUMMY CYLINDRICAL BEARING - OUTER RING 

LUBRICANT ENTERING SPHERICAL TEST BEARING 

LUBRICANT ENTERING THRUST HYDROSTATIC 
BEARING 

LUBRICANT ENTERING RADIAL HYDROSTATIC 
BEARING 

LUBRICANT IN RESERVOIR 

LUBRICANT ENTERING DUMMY CYLINDRICAL 
BEARING 


TABLE 6 

DICTIONARY OF NODES USED IN SYSTEM MODEL , EXAMPLE 1 
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ERROR MESSAGE FROM EQUATION SOLVING ROUTINE 

ALGORITHM NO. , LOOP = 

DIVERGED ITERATION - TOTAL NUMBER OF DIVERGED 
ITERATIONS = 


A) DIVERGENCE MESSAGE 

ERROR MESSAGE FROM EQUATION SOLVER 

LOOP = , ALGORITHM = , 

REQUIRED NORM'* 0 = 

CURRENT NORM = 


CURRENT NORM IS NOT WITHIN 1001 OF USER SPECIFIED 
REQUIREMENT AND NO FURTHER IMPROVEMENT CAN BE MADE, 
EXECUTION STOPS 

B) ABORT MESSAGE 


TABLE 7: TWO DIAGNOSTIC MESSAGES GENERATED 

IN THE NEWTON- RAPHSON EQUATION SOLVER 


20 


Required Norm is 


the user input value "TOL" 
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FIGURE 1 - SPHERICAL ROLLER BEARING 
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FIGURE 2: USE OF FREE FORMAT TO SPECIFY INPUT DATA. 
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NOTE: PH II REFERS TO FIRST ROLLER IN ROW #1 

PHI2 REFERS TO FIRST ROLLER IN ROW #2 
(SHOWN DOTTED) 

FIGURE 4: BEARING COORDINATE FRAME. 
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FIGURE 5: USER INPUT BEARING CLEARANCES. 
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FIGURE 9: IIIQ1 SPEED PLANET GEAR LOADING . 
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FIGURE 10: PLANET GEAR OUTER RING CROSS-SECTION 

(Effective area for computation of 
XI shown cross-hatched.) 
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FIGURE 12: TOOTH LOADS ON PLANET GEAR 
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NOTE: 



A AND B ARE INSIDE AND OUTSIDE RADII OF SPLIT INTER- 
FACE FOR A TWO PIECE CAGE. 


FIGURE 14: SPLIT CAGE INPUT DATA. 
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FIGURE lb: OPTIONAL SYMMETRIC ROLLER GEOMETRY INPIFI DATA 
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FIGURE 17: OPTIONAL ROLLER GEOMETRY DATA. 
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FIG. Al: USER INPUT DATA FOR COMPUTATION OF BEARING OPERATING CLEARANCE (OPTION 1) 
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Flfi. A3 : (IS I: K SPECIFIED SLICE WIDTH CARD FORMATS (OPTION 3) 



User may input slice widths in inches (METRIC * .FALSE, 
or millimeters (METRIC - .TRUE.) 







User may input radii in inches (METRIC » .FALSE.) or 
millimeters (METRIC = .TRUE.). 
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Hi;. AS (CONTINUED) USER SPECIFIED, COMPLETELY VARIABLE, ROLLER GEOMETRY CARD FORMATS (OPTION S) 
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OAK It 1 - Tl MI'liKAIUKC CALCULATION CONTROL CARD 
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nr.URI A6 (CONTINUED): INPUT DATA CARD FORMATS FOR TEMPERATURE CALCULATION (OPTION 6) 









HOJkfc '0 (CONTINUED): INPUT DATA CARD FORMATS TOR Tl MI'hRATlIRF. CALCULATION (OPTION 6) 
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CARD JA - NODES WHERE •EARING HEAT IS GENERATED. USE THIS CARD EUR 



HGURE A6 (CONI INUED): INPUT DATA CARD FORMATS FOR TEMPERATURE CALCULATION (OPTION 0) 



















AT81D007 




I 

I 




I 

r 

r 

t 

r. 



ORIGINAL FA'- T : F* 
OF POOR QUALITY 


\ 


102 


I'UiUKH A 6 (CONTINUED) : INPUT DATA CAKI) FORMATS FOR TEMPERATURE CALCULATION (OPTION u) 
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FIGURE A6 (CONTINUED): INPUT DATA CARD FORMATS FOR TEMPERATURE CALCULATION (OPTION 6) 















































CARD 8 - NODAL I It AT CAPACITIES (USE ONLY I'OK TRANSIENT CALCULATIONS) 


I 

( 


i. 

r 

r 



AT81D007 


ORlG?NAi- 1 j 

OF POOR QUALITY. 






AT81D007 


APPENDIX B : 

HEAT TRANSFER COMPUTATION NOTES 


SKF TECHNOLOGY SERVICES 

oKf industries. in-. 


106 


AT81D007 


B.l BASIC EQUATIONS* 


B.1.1 Heat Conduction 


The rate of heat flow q c ^ ^ (W) that is conducted from node 
i to node j may be expressed by, 


l ci , J 




(t i - V 


and 


t. are the temperatures at i and j, respectively, A. . the 

1 J 2 1 » J 

area normal to the heat flow, (m ) the distance (m) and X^ 

the thermal conductivity between i and j , (W/m°C) . 


Assuming that the structure between point i and j is com- 
posed of different materials, an equivalent heat conductivity may 
be calculated as follows: 



The calculation of the areas will be discussed in Section B.l. 5. 


*This Appendix is based on the material presented in Reference 5. 
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B.1.2 Convection 

The rate of heat flow that is transferred between a solid 
structure and air by free convection may be expressed by 


l vi,j 


a . .• 
i. j 


a. : 
i.j 


c - 


1.25 


SIGN (t. - tj) 


where 


SIGN 


1, if (t t 
-1, if (tj 


tj) > 0 

tj) < 0 


in which 

/ 2 . 5 


• 10" 2 W/m 2 - (degC)** 2 ^ for hot surfaces facing 
upward and cold surfaces facing downward 




1.4 


\1.S 


10' 2 W/m 2 - (degC) 1,25 


for hot surfaces facing 


downward and cold surfaces facing upward 
• 10’ 2 W/m 2 - (degC)^* 2 ^ for vertical surfaces 


For other special conditions, a^j must be estimated by referring 
to heat transfer literature. 

The rate of heat flow that is transferred between a solid 
structure and a fluid by forced convection may be expressed by 


q. . = ci • .A. . (t • - t • 
4 ni,j i,j i,J 1 i J 


) 


in which is the convective heat transfer coefficient. 
Now, with a * , introduce the Nusselt number 


N 


= 

U A 


the Reynolds number 

R e 

and the Prandtl number 


n _ UL 

R e v 


pv C p 

■nr 
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where 

L is a characteristic length which is equal to the diameter in 
the case of a cylindrical surface and is equal to the plate 
length in case of a flat surface (m) . 

U is a characteristic velocity which is equal to the difference 
between the fluid velocity at some distance from the surface 
and the surface velocity (m/sec) . 

A is the fluid thermal conductivity (W/M°C) 

? 

v is the fluid kinematic viscosity (M /sec) 

p is the fluid density (kg/m 3 ) 

Cp is the fluid specific heat (J/kg°C) 

For given values of R g and P r the Nusselt number N u and 
thus the heat transfer coefficient may be estimated from one of 
the following expressions: 

Laminar flow along a flat plate: R < 2300 

N u - 0. 323 £ • 

Laminar flow of a liquid in a pipe: 

N u - 1.36 V R e • P r (£) 

where D is the pipe diameter and L the pipe length 

Turbulent flow of a liquid in a pipe: 

N = 0.027 • R°' 8 • 3 ^7 
u e r 


Gas flow inside and outside a tube: 
0-3 R °- 57 


Liquid flow outside a tube: 


N - 0.6 R°' 5 
u e 


,0.31 
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Forced convection from the outer surface of a rotating 
shaft 

N * 0.11 [ 0.5 R 2 • P ] 0,35 
u e r 

where the Reynolds number R is developed by the shaft 
2 e 

R , D 

e v 


in which w is the angular velocity (rad/sec) 

D is the shaft diameter (m) 

The average coefficient of forced convection to the lubri- 
cating oil within a rolling contact bearing may be approximated 
by, 

a = 0.0986 [1 ♦ }*5x(P r ) 1/3 

m 


using + for outer ring rotation 
- for inner ring rotation 

in which w is the bearing operating speed (rad/sec) 

D is the diameter of the rolling elements (mm) 

d m is the bearing pitch diameter (mm) 

8 is the bearing contact angle; zero for cylindrical 
roller bearings (degrees) 

B.1.3 Fluid Flow 

The rate of heat flow that is transferred from fluid node 
i to fluid node j by fluid flow is 


fi, j 


pv ij C p (t. - tp 


is the volume rate of flow from i to j . It must be observed 
that the continuity of mass requires the following equation to 
be satisfied 

■ 0 

provided the fluid density is constant. The summation should be 
extended over all nodes i within the fluid which have heat 
exchange with node j by fluid flow. 
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B.1.4 Heat Radiation 

The rate of heat flow that is radiated to node j from 
node i is expressed by 

iRi.J ’ 6 iJ * 273 > 4 - * 273 ’ 4} 

where 

Tj - t. * 273.16 
T i - t i ♦ 273.16 


and the value of the coefficient 5. . depends on the geometry and 

1 » J 

the emissivity or the absorptivity of the bodies involved. 

For radiation between large, parallel and adjacent sur- 
faces of equal area, A. . and emissivity, e. . , 6. . is obtained 

i»j 1 > J 1 » J 

from the equation 




where o, the Ste fan-Boltzmann constant, is 
a - 5.76 • 10' 8 W/m Z /(degK) 4 


For radiation between concentric spheres and coaxial 
cylinders of equal emissivity, e. . , <5 ■ - is given by the 

* 9 J 1 9 J 

equation 


±L± 


a A. 


1 


ij \ , 

1 + Ci-e 

1,J A i » j 


where a is as above A. . is the area of the enclosed body and 

1 * 1 

A*. . is the area of the surrounding body, i.e., A. . < A*. .. 

i , J 1 » J 1 * J 


Expressions for <5- . that are valid for more complicated 

1 » J 

geometries or for different emissivities may be found in the 
heat transfer literature. 
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B.1.5 Calculation of Areas 

In the case of heat transfer in the axial direction A. . 

1 f J 

is given by the equation (Fig. B-3) 

A. . * 2IIr • Ar 
i » J ni 

Referring to the temperature calculation input instructions, 
card 7, but recalling L must be input in mm not m. 


L 7 = Ar = r 7 


- r. 


In the case of heat transfer in the radial direction, 

A. . is obtained from the expression 
1 » J 


i.J 


2 II r_ 


m 


H; L. 


r m ; L 2 


H 


and similarly for the radiation term above 
A*. . * 2IIr* H 

i,j m 



L 2 = 2H 

in which H is the length of the cylindrical surface; where heat 

is conducted between i and i, r is given by the same equation 

m 

as above (Fig B-4 (a)); where heat is convected between i and j, 
r is the radius of the cylindrical surface (Fig. 3-4(b)); where 

TTl 

heat is radiated between i and i , r is the radius of the enclosed 

J m 

cylindrical surface and r * the radius of the surrounding cvlin- 
* m 

drical surface (Fig. B-4(c)). 
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Fig . B - 3 


Fig. B’4(a) 


B. 2 TRANSIENT ANALYSIS 

For the transient analysis all of the data pertaining to 
the node to node heat transfer coefficients must be provided by 
the input. Additionally, the volume and the specific heat at 
each node is required. For metal nodes this input is straight- 
forward. To date, fluid nodal volume in a r ree space such as 
the bearing cavity, has not been correlated with lubricant flow 
rate or application method. 


Fig. B -4 (b) 


Fig. B *4 (c) 
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PROGRAM OUTPUT, BXAMPI.F 1: PAGE 14 
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FROCK \M OUTPUT, EXAMPLE 1: PACE 
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PROGRAM OUTPUT, EXAMPLE 1: PAGE 18 
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PROGRAM OUTPUT, PXAMPI F 1: PAGE 19 
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PROGRAM OUTPUT, EXAMPLE 1: PAGE 20 
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PROGRAM GU I PUT, EXAMPLE 1: PAGE 
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PROGRAM OUTPUT, EXAMPLE 1: PAGE 22 
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PROGRAM OUTPUT, EXAMPLE 1: PAGE 23 
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PR0f?R\M OUTPUT, EXAMPLE 1: PAGE 24 
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PROGRAM OUTPUT, EXAMPLE 2: PAGE 
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PROGRAM OUTPUT, EXAMPLE 2: PAGE 



